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ABSTRACT 

The recent discovery of burst oscillations at 1122 Hz in the x-ray transient XTE J1739-285, together 
with the measurement of the mass of the binary millisecond pulsar PSR J0751+1807 (2.1 ± 0.2 M Q ) 
can finally allow us to put strong, model-independent observational constraints to the equation of state 
of compact stars. We show that the measurement of the moment of inertia of PSR J0737+3039A, 
together with these constraints, could allow to discriminate further the details of the inner structure 
of neutron stars. Moreover, we show that if XTE J1739-285 is constituted of nucleonic matter, any 
equation of state allows only a narrow range of very high masses, and this could explain why up to 
now compact stars spinning faster than a millisecond have been so difficult to detect. 
Subject headings: Stars: neutron - X-rays: binaries - binaries: close - relativity. 



1. INTRODUCTION 

Ever since the discov ery of the first neutron star (NS) 
bv lHewish et alj (|1968f ). these compact objects have been 
regarded as the ideal experimental test-bed for the the- 
ories of the state of matter at supernuclear densities. In 
particular, the mass, the radius, the attainable spin fre- 
quency of a NS strongly depend on the equation of state 
(EOS) governing matter at these densities. 

Different models for the EOS of dense hadronic matter 
predict a NS maximum mass (M max ) in the range of 1.4- 
2.6 Mq, and a corresponding central density n c in range 
of 4-10 times the saturation density ( rip = 0.16 fm ) 
of nuclear matter (|Prakash et al.lfl997h . In the case of a 
star with M ~ 1.4 M Q , different EOS models predict a 
radius in the range of 7-16 km. 

Due to the large value of the stellar central densities, 
various regimes (particle species and phases) of dense 
hadronic matter are expected in the interiors of NS. Con- 
sequently, different types of "neutron stars" (compact 
stars, CSs) are hypothesized. In the simplest model the 
core of a neutron star is described as a uniform fluid of 
neutron rich nuclear matter in equilibrium with respect 
to the weak interaction: these are the so-called "tradi- 
tional" NSs. However, due to the large value of the stellar 
central density and to the rapid increase of the nucleon 
chemical potentials with density, hyperons (A, E _ , £ , 
E + , S~ and S° particles) are expected to appear in the 
inner core of the star. Other exotic phases of hadronic 
matter such as a Bose-Einstein condensate of negative 
pion (ir~) or negative kaon (K~) could be present in the 
inner part of the star. The core of the more massive NS 
is also one of the best candidates in the Universe for a 
phase transition from hadronic matter to a deconfincd 
quark phase to occur. Compact stars which possess a 
"quark matter core" either as a mixed phase of decon- 
fined quarks and hadrons or as a pur e quark matter (QM ) 
phase are called Hybrid Stars (HyS, Glcndcn nindll996] ). 
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The more conventional NSs in which no fraction of QM 
is present, are referred to as pure Hadronic Stars (HSs). 

Even more challenging than the existence of a quark 
core in a NS, is the possible existence of a new type of 
CS consisting completely of a charge neutral deconfincd 
mixture of up, down, strange quarks and electr ons, sat- 
isfying the hypoth esis on the absolute stability (jBodmerl 
U97lMWitten1ll984D of strange quark matter (SQM). Such 



CSs have been called strange stars (SS). The analysis of 
different type of observational data has given in direct ev- 
idenc e for the possible existence of SS (see e.g. (|Li et alJ 
1999)). We will refer to hybrid stars and strange stars 
collectively as Quark Stars (QS). 

Unluckily, for decades no significant constraint on both 
the mass and the radius of a single object could be ex- 
tracted from astrophysical observations. While the de- 
termination of the masses of some CSs was made pos- 
sible thanks to the study of the timing motion of ra- 
dio pulsars in binary systems, the radius determina- 
tion, or alternatively the moment of inertia, has been 
more elusive. Today the measured masses of CSs lie 
around the average value of 1.35 ± 0.04 M Q and are 
roughly distributed according t o a Gaussian distribution 
(jThorsett fc Chakrabartvlll999l ). 

In recent years the discovery of kilohertz quasi-periodic 
oscillations (QPOs) in the power spectra of several ac- 
creting CSs has opened new possibilities, as they are 
phenomena whose timescale corresponds to the dynam- 
ical timescale in the very neighborhood of a CSs (see 
Ivan der Kl is 2006, for a review). These QPOs have peak 
frequencies in the 200-1200 Hz range and usually appear 
as two broad peaks (denoted upper and lower kHz QPO), 
that simultaneously move up and down according to the 
source accretion state. The lower kHz QPO sometimes 
becomes undetected as the upper kHz QPO reaches lower 
frequencies; for low accretion rates both peaks can disap- 
pear. While their peak frequencies smoothly change with 
the accretion rate of the source, their peak separation is 
roughly constant and equal, within 20%, to the spin fre- 
quency of the NS, or for some sources, to the half of it 
(|Wijnands et al.ll2003h . The reason for this discrepancy 
is still unaccounted for. Most theoretical models of kHz 
QPOs associate the frequency of the highest peak to the 
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motion of matter at the last stable orbit of an accretion 
disk around the CS. If this is true, assuming that matter 
in the accretion disc is in Keplerian motion, it is possible 
to derive an upper limit to the density of the CS. This, 
in tur n, puts some limits to the EOS of CSs (see for ex- 
ample Zh ang et al.ll2006T ). but still relies on a theoretical 
interpretation of a phenomenon which is not completely 
understood. 

Another observational constraint comes from the re- 
ported gravitational redshift of absorption lines produced 
in th e photosphere of the X-ray binary system EXO 0748- 
676 (|Cottam. Paerels fc Mendel l2002t ) . The redshift (z 
= 0.35) implies a compactness parameter for the CS of 
GM/Rc 2 ~ 0.23. However, even this result (which any- 
way is not able to constrain the EOS significantly if taken 
alone) is dependent both on the model adopted for the 
spectral fitting and on the assumptions on the emitting 
region. 

Another astrophysical observable is the radiation ra- 
dius Rao of the emission from isolated CSs: the nearest 
and better studied isolated CS, RX J1856. 5-3754, has 
a spectrum that can be well fitted with a blackbody 
of 57 eV temperature. This value sets an upper limit 
for the Rgo of this object, which was estimated < 8 km 
(|Drake et al J 12002). This limit, which is quite demand- 
ing, can imply an EOS based on QM. However, the de- 
termination of this radius is strongly dependent on the 
distance of the source and on theoretical issues regarding 
the geometry of the emission and t he composition of the 
atmosphere. iTrumper et al.l (|2004l ) has recently argue d 
against the constraints found by iDrake et al.l {2002), 
claiming that the upper limit is much higher. 

The situation is even worse if we try to infer the ra- 
dius of the CS from the blackb ody temperature o f an 
accreting X-ray binary (see e.g. IZhang et al.l l2006h : in 
this case even the origin of the emission can be debat- 
able, as the blackbody emission can be originated near 
the stellar surface as well as from the accretion disc. 

A truly model-independent constraint has been found 
with the measurement of the mass of several millisec- 
ond radio pulsars thanks to the detection of post- 
New tonian corrections to their orbital motion. In partic- 
ular, iNice et alJ (|2005f ) found that PSR J 0751+1807 has 
a mas s of 2.1 ± 0.2M Q (1 a error), and iRansom et aTl 
( 2005) found that pulsar I of the globular cluster Terzan 
5 has a mass that exceeds 1.68 iW© at 95% confidence 
level, showing that CSs can have quite large masses. 

Another model independent limit to the EOS at su- 
pernuclear densities comes fro m laboratory experi ments 
on heavy ions collisions (HIC): iKlahn et all (|2006l ) show 
that some constraints on the EOS can be given by flow 
data and sub-threshold Kaon production in HIC. They 
found that these c onstraints would rule out some of the 
very stiff EOSs. iGrigorian. Blaschke fc Klahnl (|2006l ) 
solved the problem by simply assuming that above a cer- 
tain energy density, a quark core forms, with a phase 
transition to stiff, color superconducting QM that can 
alter only the behavior of the EOS at very high densities, 
making it softer near the maximum mass. In practice, 
this results in small modifications to the behavior of the 
EOS that do not alter very significantly the mass-radius 
relationship of CSs. 

Finally, we will show that the recent discovery of burst 
oscillations with a frequency of 1122 Hz in XTE J1739- 



285 (|Kaaret et al.l |2006[ ). is potentially another im- 
portant constraint to the EOS of CSs: these oscil- 
lations, occurring during type I X-ray bursts, have 
the same frequency of t he spin frequency of the CS 
(|Chakrabarty et all l2003h . This discovery can put a 
real, model independent limit to the radius of the NS 
in XTE J1739-285, and thus to its EOS. 

This discovery rules out some speculations that were 
made in order to explain the inability of CSs to spin 
up to periods below one millisecond. It was argued 
that the lack of these fast spinning objects was due 
to t he emission of gravitational waves during accre- 
tion (lAndersson. Kokkotas fc Stergioulaslfl998t iBildstenl 
Il998h or to the disc- ma gnetic field interaction in low mass 
X-ray binaries (see e.g lAndersson et al.l l2005). All these 
theories have been wiped out by this discovery, but the 
problem has now turned to: why only one CS is observed 
to spin below one millisecond? 

In this letter we first impose the model independent 
constraints to the EOS of CSs, we discuss how future 
further observations can improve these constraints, then 
we try to explain the difficulty in finding very fast objects 
on the basis of our findings on the structure of CSs. 

2. MODEL-INDEPENDENT CONSTRAINTS 

In figure [T] we show in the mass-radius plane the se- 
quences for non-rotating CSs with very different com- 
positions, together with constraints coming from model- 
independent observational results. First, the sequence 
should include masses as large as the one of PSR 
J0751+1807 (horizontal dashed lines) and of pulsar I of 
Terzan 5 (dotted line). We assume that the maximum 
non-rotating mass is almost equal to the maximum mass 
at the spin frequencies of these objects, which is a good 
approximation given the relatively low spin frequency of 
both objects. Second, the mass-radius relation should be 
such to allow the star to spin up to the maximum spin 
rate observed to date, i.e. to 1122 Hz, as observed in 
XTE J1739-285. 

The maximum rate of rotation (f2 m ax) sustain- 
able by a CS strongly depends on the overall 
stiffness of the EOS of dense hadronic matter. 
Equilibrium sequences of rapidly rotating CSs have 
been constructed numeri cally in general relativ- 
ity by sev eral groups (ICook. Shapiro fc Teukolskyl 

1994; iDatta. Thampan fc Bombacil I199S : 

Bombaci. Thampan fc Dattal |200*oT k The numerical 
results for f2 max obtained for a broad set of realistic 
EOS can be reproduced with a very good accuracy using 
simple empirical formulas which relate f2 ma x to the mass 
(M m „) and radius Rp of the non-rot ating maximum 
mass configuration ([Salgado et al.|[l994T ). The following 
simple empirical formula given by Lattimer (Conference 
"Isolated Neutron Stars, London 2006, see also Latimer 
& Prakash 2004) approximately describes the minimum 
rotation period P m i n = 27r/f2 max for a star of mass M 
and non-rotating radius R 
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Using the rotation period P = 0.891 ms, correspond- 
ing to the measured X-ray burst oscillation frequency for 
XTE J1739-285, we get the following upper limit for the 
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Fig. 1. — Mass-radius plane with mass limits (dashed black lines) 
and spin and causality limits (green lines), plus sequences for all 
EOSs considered (see text for details). 

radius of the compact star in XTE J1739-285 

R < 9.52 (M/M Q ) 1/3 km. (2) 

This condition is plotted in figure as a green line, with 
another green line indicating the causality limit for the 
mass-radius relation. 

As we said before, there exist various possible classes 
of CSs. We considered several EOSs in order to match 
them with observations, and constructed non-rotating se- 
quences varying the central density of the star. In fig- 
ure [T] we report the mass-radius sequences for several 
EOSs. The curves relative to NSs are colored in black, 
those for hyperon stars in brown. In blue we show the 
sequences for HySs, and in red those for SSs. In par- 
ticular, we report results for the following classes of sys- 
tems: i) pure neutron matter CSs (L) calculated within 
a relativistic field theoretical approach in the mean field 
appro ximation with meson exchange (|Arnett fe Bowersi 
I1977J ). This outdated and schematic model, is consid- 
ered as an example of an extemely stiff EOS; ii) nucle- 
onic stars (label BBB2) built with an EOS calculated 
using the microscopic Brueckner-Hartree-Fock (BHF) 
many-body approach with the Paris nucleo n-nucleon in- 
terac tion plus nuclear three-body forces (jBaldo et al.1 
119971) ; iii) nucleonic stars (KS) bu ilt with an EOS de- 
rived (jKrastev fc Sammaruc a 2006) from the Dirac-BHF 
approach using the Bonn B nuclear interaction, iv) 
nucleonic stars (GM1) calculated within a relativistic 
field theoretical approach in the mean field approxima- 
tion (|Glendenning fc Moszkowk"llll991l ): v) hyperon stars 
(Hyp) calculated with the previous GM1 EOS with the 
inclusion of hyperons; vi) Hybrid star sequences cal- 
culated using the same GM1 EOS for hyperonic mat- 
ter to descri be the hadron i c pha se, and the MIT bag 
model EOS (jFarhi fc Jaffd 1 1984 ) to describe the de- 
confincd quark phase. For the masses of the three 
quark flavors we took m u = rrid = and m s — 150 
MeV. For the value of the bag constant we took B = 
208.24 MeV/fm 3 (Hybl) and B = 136.63 MeV/fm 3 
(Hyb2). The curve (Hyb3) refers to hybrid stars con- 
str ucted using a different parametriza tion for the EOS 
by iGlendenning fc Moszkowkl (fl99lt ) (GM3) for the 
hadronic phase and u sing B = 80 MeV /fm 3 ; vii) Hy- 
brid stars (Hyb(CFL), lAlford et aT1l2005h whose core is a 
mixed phase of nuclear matter and color-superconducting 
quark phase in the so called "color flavor locked" (CFL) 



phase; viii) strange star sequences (B85, B70) built with 
the MIT bag model with with B = 85 MeV/fm 3 and B = 
70 MeV/fm 3 respectively {m u = m,d — and m s = 150 
MeV), and sequence (B60) with B = 60 MeV/fm 3 and 
massless quarks; ix) CFL s trange stars (label SS(CFL)) 
(jLugones fc Horvathl [2003] ) within the bag model EOS 
with B = 70 MeV/fm 3 , m u = m d = and m s = 150 
MeV and quark pairing gap A = 100 MeV. 

Looking at the figure we see that, although almost all 
classes of EOSs are still allowed (apart probably from 
Hyperon stars), only the stiffest SSs are allowed, and that 
while soft HSs are excluded by the mass constraint, some 
of the stiffest ones (i.e. Hybl and Hyb2, L) are excluded 
by the rotational limit. Other EOSs are only marginally 
compatible with the constraints (BBB2, GM1, B60 and 
CFL) and would require exceptional conditions (e.g. a 
mass almost equal to the limiting mass for gravitational 
collapse) to meet the observational constraints. 

3. DISCUSSION AND CONCLUSIONS 

The exceptional discovery of a CS spinning with a 
period below one millisecond allowed us to put a new 
model-independent constraint to the EOS of matter 
above nuclear density. Combining this with the con- 
straint given by the highest masses measured in millisec- 
ond radio pulsars, we are able to narrow the classes of 
EOSs that could describe compact objects. As it is clear 
from figure Q] we are left with two main classes of EOSs 
that comply with our constraints: very stiff quark matter 
EOSs and some quite stiff nucleonic matter EOSs. These 
two classes have completely different characteristics. 

SSs that satisfy both constraints have radii smaller 
than ~ 10.5 km for a mass of 1.337 M (the mass of 
the ultrare lativistic pulsar P S R J7 37-303A). Using the 
formula bv lBeiger fc Haensell (|2002D . the moment of in- 
ertia of a SS with a mass of 1.4 M Q and a radius of 10.5 
km is, being x = M/R (M©/km) 

I = -(1 + x)MR 2 = 1.32 x 10 45 g cm 2 . (3) 
5 

On the other hand, HSs that comply to both con- 
straints have radii in the range 12 — 13.5 km for a mass 
of 1.4 M Q , and cons equently larger moment of inertia: 
using the formula by iBeiger fc Haensell (|2002D for HSs, 
assuming R ~ 13 km we obtain 

I =-(! + 5x)MR 2 = 1.52 x 10 45 g cm 2 . (4) 
9 

There is a 13% discrepancy between the two values, so 
that the measurement of the moment of inertia of PSR 
J0737-3039A with the uncertainty of 10% - quoted as 
reachable bv lLattimer fc Schutzl (|2005l ) - could tell us if 
this pulsar is a SS or a very stiff HS, although it will not 
be possible to discri minate between a NS and a HyS, as 
Klah n~et al.l (|2006bt ) have shown. We have to note that if 
the HS is one of the less stiff ones R ~ 11 km, its moment 
of inertia would be virtually indistinguishable from the 
one of a stiff SS. However, soft NSs (like EOS BBB2 of 
figure 1) are strongly disfavored by the mass constraint. 

If we look at figure[TJ we can also notice that the mass- 
radii sequences of the allowed HSs cross the spin limit line 
for XTE J1739-285 only at very high masses: as an ex- 
ample, the KS EOS allows only stars exceeding 2 Mq to 
attain such high spins. This is because the radii of stiff 
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HSs shrink only when approaching the maximum allowed 
masses. So, not only the moment of inertia of the HSs is 
large, meaning that it is necessary to accrete more mass 
to get a star to periods in the order of one millisecond, 
but also that only a very massive star would be able to 
reach a period of 0.891 ms, and only at mass shedding. 
Moreover, the HS should lie in a very narrow range of 
masses in order to be stable against gravitational col- 
lapse: if the star is supramassive (which can be the case 
of XTE J 1739-285), accretion of further matter can eas- 
ily le ad to the collapse to a black hole ([Lavagetto et al.l 
l200l . 

On the contrary, since SSs are self-bound objects, their 
radii expand with growing mass, allowing very fast spins 
even for lower masses (i.e. EOS B60 can attain this pe- 
riod for M > 1.5Af ). This means that it is much easier 
to attain the period of XTE J1739-285 in this case. 

Paradoxically, if - as the discovery of XTE J 1739- 
285 implies - there is no mechanis m avoiding the spin 
up to millisecond periods (see e.g. Chakra bartv et aTl 
l2003f ). SSs spinning at sub-millisecond periods should be 
quite common, as accreting ~ 0.2 M W could spin the S S 
to periods below one millisecond (IZdunik et al.l [2 002). 
Even evolutionary arguments proposed bv lBurderi et all 
(2001), that can explain the non-detection of ultrafast 
radio pulsars do not apply to the detection of burst os- 
cillations, whi ch is thought to be observationally unbi- 
ased (see e.g. Chakrab artv! [2005T ). Thus the observa- 
tional shortage of submillisecond CSs (only one out of 
19 X-ray binaries that show coherent pulsations in their 
lightcurve) would be somewhat puzzling. 

On the other hand if CSs are composed of nucleonic 
matter, the high stiffness required to the EOS by the 
high measured masses can alone explain the observa- 
tional shortage of very fast CSs: a stiff EOS implies 



that a submillisecond HS should both be very massive 
and reach the mass-shedding regime, which seem diffi- 
cult given only a few measured NS masses are larger than 
2 Mq (see e.g. lLattimer fc P rakash 2004). Moreover, as 
the range of masses that allow a stable submillisecond HS 
is very narrow (i.e. in the case of EOS KS the allowed 
mass range is 2— 2.2 Mq, and it is even narrower for other 
EOSs, sec figure!]}, this shortens further the chances we 
have to observe a submillisecond HS, as a supramassive 
HS can either collapse during accretion, or collapse due 
to braking if accretion is stopped when the star is in the 
supramassive regime: in general, a sup ramassive CS is 
expec ted to have a very short lifetime (|Lavagetto et al.l 
l2004h . We have therefore a strong selection effect against 
the detection of very fast spinning objects. Again, if 
there is a transition to deconfined quark matter at very 
high densities in NSs, as required by the results of HIC 
experiments (jKlahn et al.l l2006t iGrigorian et~ai] l2006ft . 
this will not change significantly this scenario. 

In conclusion, we showed that when the moment of in- 
ertia of PSR J0737-3039A would be measured with an 
accuracy of ~ 10%, we will likely be able, by combining 
this measurement with the constraints presented here, to 
assess the nature of CSs using only model-independent 
observational constraints. In the meanwhile, if no other 
object spinning below one millisecond is found, this is 
probably an argument favoring the hypothesis that CSs 
are composed of nucleonic matter. Anyhow, this pos- 
sibility challenges our present understanding of particle 
physics at supernuclear densities. In fact, one should 
envisage a plausible physical mechanism to avoid the ap- 
pearance of non-nucleonic degrees of freedom (hyperons, 
quark matter, etc.) in the dense interiors of neutron 
stars. 
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